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Abstract 

We analyze collocation methods for nonlinear homogeneous Volterra-Hammerstein 
integral equations with non-Lipschitz nonlinearity. We present different kinds of exis- 
tence and uniqueness of nontrivial collocation solutions and we give conditions for such 
existence and uniqueness in some cases. Finally we illustrate these methods with an 
example of a collocation problem, and we give some examples of collocation problems 
that do not fit in the cases studied previously. 

1 Introduction 

The aim of this paper is the numerical analysis of the nonlinear homogeneous Volterra- 
Hammerstein integral equation (HVHIE) 

y{t)^{Hy){t):^ f Kit,s)G{y{s))ds, teI:=[0,T], (1) 
Jq 

by means of collocation methods on spaces of local polynomials. This equation has multiple 
applications in analysis and physics, as for example, the study of viscoelastic materials, the 
renewal equation, seismic response, transverse oscillations or flows of heat (see [Illl]). 

Functions K and G are called kernel and nonlinearity respectively, and we assume that 
the following general conditions are always held, even if they are not explicitly mentioned. 

• Over K. The kernel if : — > [0, +oo[ is a locally bounded function and its support 
is in {(t, s) g . Q < s < i}. 

For every t > 0, the map s K {t, s) is locally integrable, and K {t, s) ds is a strictly 
increasing function. 

• Over G. The nonlinearity G : [0,-|-c»[ — ^ [0, +oo[ is a continuous, strictly increasing 
function, and G (0) — 0. 

*Work supported by project MTM2008-05460, Spain. 
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Note that since G{0) — 0, the zero function is a solution of equation ([T]), known as 
trivial solution, and so, uniqueness of solutions is no longer a desired property for equation 
([T]) because we are obviously interested in nontrivial solutions. Existence and uniqueness of 
nontrivial solutions of equation (IT]), as well as their properties, have been deeply studied in a 
wide range of cases for K and G [olIlllSJinilZllEli especially in the case of convolution equations, 
i.e. K{t,s) = k{t — s). In general, a necessary and sufhcient condition for the existence of 
a nontrivial solution is the existence of a nontrivial subsolution; that is, a positive function 
u such that u{t) < (%u) (t). So most of the results on existence of nontrivial solutions are, 
indeed, characterizations of the existence of subsolutions. For instance, in |5] it can be found 
the next result: under the general conditions, equation ([T]) has a nontrivial solution if and 
only if there is a positive integrable function f{x) such that JC {J^{x) — J^{s)) ds > G^^{x), 
X > 0, where IC{x) :— k{s) ds and J-{x) := f{s) ds. 

It is important to note that usually, in the analysis of solutions for non-homogeneous 
Volterra integral equations (and their numerical approximations), most of the existence and 
uniqueness theorems require that a Lipschitz condition is held by the nonlinearity (with 
some exceptions, for instance [9]). This is not our case, since it is well known that if the 
nonlinearity is Lipschitz-continuous, then the unique solution of ([T]) is the trivial one [TU] . 
Thus, the case we are going to consider in this paper is beyond the scope of classical results 
of numerical analysis of non-homogeneous Volterra integral equations, in the sense that we 
need a non-Lipschitz nonlinearity. 

Actually, there is a wide range of numerical methods available for solving integral equa- 
tions (see [llj for a comprehensive survey on the subject): iterative methods, wavelet methods 
[12l [131 [Mj [15] , generalized Runge-Kutta methods [El [17], or even Monte Carlo methods [18]. 
Collocation methods flO', have proved to be very suitable for a wide range of equations, 
because of their accuracy, stability and rapid convergence. In this work we use collocation 
methods to solve the nonlinear HVHIE ^ written in its implicitly linear form (see below). 
We also give conditions for different kinds of existence and uniqueness of nontrivial collocation 
solutions for the corresponding collocation equations. 

We organize this paper into four sections. In Section [2] we write equation ([I]) in its 
implicitly linear form and we describe the corresponding collocation equations; moreover, we 
define the concept of nontrivial collocation solution. In Section [3| we present different kinds 
of existence of nontrivial collocation solutions and we give conditions for their existence and 
uniqueness in some cases, considering convolution and nonconvolution kernels. In Section 
|4| we illustrate the collocation methods and their numerical convergence with an example, 
showing how the errors change as the collocation points vary. Moreover, we give some 
examples of collocation problems that do not fit in the cases studied in the paper. Finally, 
we present the proofs of the main results in an appendix at the end of the paper, for the sake 
of readability. 



2 Preliminary concepts 

Let us consider the nonlinear homogeneous Volterra-Hammcrstein integral equation (HVHIE) 
given by ([I]). Taking z :— G o y, equation (jl]) can be written as an implicitly linear homoge- 
neous Volterra integral equation (HVIE) for z: 

z{t)^G{iVz){t)) = G(^J^ K{t,s)z{s)dsy tel, (2) 

where V is the linear Volterra operator. So, if z is a solution of (|2|, then y :— Vz is a solution 
of Q. It is known (see [IHl p. 143]) that, under suitable assumptions on the nonlinearity G, 
there is a one-to-one correspondence between solutions of ([T]) and ([2|. Particularly, if G is 
injective then y = G~^ o z and hence this correspondence is given, which is granted by the 
general conditions exposed above. 



2 



2.1 Collocation problems for implicitly linear HVIEs 

First, we are going to introduce the collocation problem associated to equation ([2]), and give 
the equations for determining a collocation solution, that we will use for approximating a 
solution of ([2| or ([l]) (see Remark 2.1 below). 

Let Ih '■— {tn '■ — to < ti < . . . < tjf ~ T} be a mesh (not necessarily uniform) on the 
interval / = [0, T], and set (t„ :— i^+i] with lengths /i„ := — tn {n = 0, . . . , N — 1). 
The quantity h :— max{/i„ : < n < A'^ — 1} is called the stepsize. 

Given a set of m collocation parameters {ci : < ci < . . . < c,„ < 1}, the collocation 
points are given by t^ i := tn + Cih^ {n — 0, . . . , — 1) (i = 1, . . . , m), and the set of 
collocation points is denoted by X^. 

All this defines a collocation problem for equation ^ (see [3T], (TUl p. 117]), and a 
collocation solution Zh is given by the collocation equation 

Zh (t) ^g(^J^K {t, s) Zh (s) ds^ , te Xh , (3) 

where Zh is in the space of piecewise polynomials of degree less than m (see [101 p. 85]). Note 
that the identically zero function is always a collocation solution, since G (0) = 0. 

Remark 2.1 From now on, a "collocation problem" or a "collocation solution" will be always 
referred to the implicitly linear equation ([2]). So, if we want to obtain an estimation of a 
solution of the nonlinear HVHIE ([I]), then we have to consider yh := Vzh- 

As it is stated in [TU] , a collocation solution z^ is completely determined by the coefficients 
ZnA Zh {tn.i) (n = 0, . . . , iV - 1) (i = 1, . . . , m), since Zh {tn + vK) = (") ^"^J 

for all V E ]0, 1], where Lj (v) :— YiT^j (j = 1; • • ■ ^re the Lagrange fundamental 

polynomials with respect to the collocation parameters. The values of Zn^i are given by the 
systems 

Zn,^ = G {^n {tn..i) + hnj^ ^» ^^^^ ' C^) 

where ^ 

Bn{i,j).= / i4r(t„,i,t„ + s/i„)Lj (s) ds. (5) 



and ^ 

F„ (t) := r K {t, s) Zh is) ds. (6) 
Jo 

The term (^n,i) is the lag term, and can be expressed in the form 

n— 1 m 
1=0 i=l 

where ^ 

Bl{i,j):= / is:(t„.„<, + s/t()ij (s) ds 
Jo 

with n = 0, . . . , N — 1, / = 0, — 1, i = l,...,m, j — 1,...,to. 

Remark 2.2 For convolution kernels, K {t, s) ~ k{t — s), expression ([s]) is given by 

Bn{i,j)= / A: ((ci - s) /i„)ij (s) ds. 



In this case, i?„ (i,i) is independent from i„ and, given some collocation parameters, it only 
depends on hn- 
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Remark 2.3 The coefRcients Zn,i > given by Q are positive, since G is a positive function. 
But it docs not imply that Zh was positive. 

The advantage of imphcitly hnear collocation methods (called new collocation-type meth- 
ods, see [lOl p. 118]) lies in the fact that, in contrast to direct collocation methods for ([T]), 
the integrals need not to be re-computed for every iteration step when solving the nonlinear 
algebraic system Q. 



2.2 Nontrivial collocation solutions 

In this section we are going to recall the definition of nontrivial solution for the implicitly 
linear HVIE ([2| with convolution kernel, and its corresponding collocation problem. Nev- 
ertheless, it can be easily extended for the original problem given by equation ([ij, and for 



nonconvolution kernels (see Remark 2.4 below), but first we need the following definition: 



Definition 2.1 We say that a property V holds near zero if there exists e > such that V 
holds on ]0, 5[ for all < 5 < e. On the other hand, we say that V holds away from zero if 
T holds on ]i, -foo[ for ah t > 0. 

Given an implicitly linear HVIE ([2| , the zero function is always a solution, as it happens 
with equation Q. Moreover, for convolution kernels, given a solution z{t) of ([2|, and < 

c < T the c-translated function of z given by 



Zc (t) 



if < t < c 

z{t~c) if c<t<T 



is also a solution of ([2]). Thus, for convolution kernels, we say that a solution is nontrivial if 
it is neither identically zero nor a c-translated function of another solution. 

In this case, z is nontrivial if and only if it is not identically zero near zero. This charac- 
terization allows us to extend the concept of nontrivial solution to collocation problems with 
convolution kernels: 



Definition 2.2 Given a collocation problem with convolution kernel, we say that a colloca- 
tion solution is nontrivial if it is not identically zero in (Jq. 

Remark 2.4 The concept of nontrivial collocation solution can be easily extended to non- 
convolution kernels. Nevertheless, we have to take into account that the c-translation of a 
solution of an implicitly linear HVIE with nonconvolution kernel is not necessarily a solution, 
and there can exist solutions that are c-translations of functions which are not solutions. 

3 Existence and uniqueness of nontrivial collocation so- 
lutions 

Given a kernel K, a nonlinearity G and some collocation parameters {ci, . . . , Cm}, our aim 
is to study the existence of nontrivial collocation solutions (of the corresponding collocation 
problem) in an interval / = [0,T] using a mesh I^. We are only interested in existence (and 
uniqueness) properties where h can be arbitrarily small and N arbitrarily large, because 
collocation solutions should converge to solutions of ^ (if they exist) when /i — > 0+ (and 
— ^ -l-oo). Unfortunately, since uniqueness is not guaranteed, these convergence problems 
are, in general, very complex and they are not in the scope of this work. Taking this into 
account, we are going to define three different kinds of existence of nontrivial collocation 
solutions. 
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• We say that there is existence near zero if tlicre exists Hq > such that if < 
ho < Hq then there are nontrivial collocation solutions in [0,ti]; moreover, there exists 
Hn > such that if < hn < Hn then there are nontrivial collocation solutions in 
[0, tn+i] (for n = 1, . . . , A'^ — 1 and given ho, ... , hn-i > such that there are nontrivial 
collocation solutions in [0,i„]). Note that, in general, depends on /ig, • ■ ■ 

This is the most general case of existence that we are going to consider. Loosely 
speaking this definition means that "collocation solutions can always be extended a bit 
more". Nevertheless, the existence of nontrivial collocation solutions for arbitrarily 
large T is not assured (for instance, if there is a blow-up). So, this is an existence near 
zero. 

• We say that there is existence for fine meshes if there exists H > such that 
a < h < H then the corresponding collocation problem has nontrivial collocation 
solutions. 

This is a particular case of existence near zero, but in this case it is assured the existence 
of nontrivial collocation solutions in any interval /, using fine enough meshes Ih- 

• We say that there is unconditional existence if there exist nontrivial collocation 
solutions in any interval / and for any mesh Ih. 

We are going to study two cases of collocation problems: 

• Case 1: m = 1 with ci > 0. 

• Case 2: m — 2 with ci = 0. 

In these cases, the system Q is reduced to a single nonlinear equation, whose solution is 
given by the fixed points of G (a + /Sy) for some a, /3. 

Let us state some lemmas which will be needed. Taking into account the general conditions 
over G, these results can be easily proved (see Appendix A). 

Lemma 3.1 The following statements are equivalent to the statement that is unbounded 
(in ]0, +oo[j; 

(i) There exists I3q > such that G {l3y) has nonzero fixed points for all < /3 < Pq. 

(ii) Given A> 0, there exists Pa > such that G (a + Py) has nonzero fixed points for all 
< a < A and for alio < P < Pa- 

Lemma 3.2 // is unbounded near zero but it is bounded away from zero, then there exists 
/?□ > such that G (a + Py) has nonzero fixed points for all a >0 and for all < P < Pq. 

If, in addition, is a strictly decreasing function, then there exists a unique nonzero 

fixed point. 

Lemma 3.3 // is unbounded near zero and there exists a sequence {yn\^i'^i of posi- 
tive real numbers and divergent to +oo such that lim„_j.+oo = 0, then G (a + Py) has 
nonzero fixed points for all a > and for all p > 0. 

If, in addition, ^^^^ is a strictly decreasing function, then there exists a unique nonzero 
fixed point. 

Remark 3.1 The nonzero fixed points are always strictly positive, since G is strictly positive 
in ]0, +cxd[. If, in addition, G (a + Py) has fixed points for a > 0, then these fixed points are 
necessarily nonzero and, hence, strictly positive. 

So, we are going to study different kinds of existence and uniqueness of nontrivial collo- 
cation solutions by means of the equations Q considering two special cases. 
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3.1 Case 1: m = 1 with ci > 

First, we shall consider m — 1 with ci > 0. Note that if m = 1 with ci = 0, then the unique 
collocation solution is the trivial one. 
We have the equations 

Zn.i=G{Fn{tn,i) + hnBA,i) {n = 0, . . . , N - 1) , (7) 

where ^ 

B„ :=B„(1,1)= / K{t„^i,tn + shn) ds, (8) 
Jo 

and the lag terms Fn {tn,i) are given by (|6| with i — 1. From the general eonditions imposed 
over K it is assured that i3„ > 0. 

Remark 3.2 Since K is locally bounded, we have — >■ when — > 0+. 

Now, we are in position to give a characterization of the existence near zero of nontrivial 
collocation solutions. 

Proposition 3.1 Let K be a kernel such that K{t,s) < K {t' , s) for all < s < t < t' . 
Then there is existence near zero if and only if ^^^^ is unbounded. 

Next, we are going to give some sufficient conditions for the existence and uniqueness of 
nontrivial collocation solutions. 

Proposition 3.2 // is unbounded near zero but it is bounded away from zero, then there 
is existence near zero. 

If, in addition, ^^^^ is a strictly decreasing function, then there is at most one nontrivial 
collocation solution. 

For convolution kernels K(t, s) = k{t — s), we can assure existence for fine meshes. 



Proposition 3.3 // is unbounded near zero and there exists a sequence {yn}nTi 

positive real numbers and divergent to +oo such that lim„_j.+oo '^^y"^ — 0, then there is 
unconditional existence. 

If, in addition, is a strictly decreasing function, then there is at most one nontrivial 
collocation solution. 



In the proofs of Propositions 3.1 |3.2| and |3.3| (see Appendix A) it is shown that the 
nontrivial collocation solutions Zh are strictly positive. Moreover, if K{t,s) < K{t',s) for 
all < s < i < t' and taking into account ([3|, it can be easily proved that these collocation 
solutions are strictly increasing functions. 

3.2 Case 2: m = 2 with ci = 

Now, we are going to consider m = 2 with ci = 0. Hence, we have to solve the following 
equations: 

ZnA = G(F„(<„,i)) (9) 
Zn,2 = G(F„(t„,2) + /inBn(2,l)Z„,i + /i„B„(2,2)Z„,2) (n = 0, . . . , TV - 1) , (10) 

where 

rc2 



B„(2,j)== / K{t„,2,tn + sh„)Lj{s) ds, (j = l,2) (11) 
Jo 

and Fn {tn,i) is given by ([6) with m — 2. Note that for solving Zn,2 this system of equations 
can be reduced to the single equation (10), since ([9| gives us Zn,i directly. 
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From the general conditions imposed over K and taking into account that functions 
Li (s) = 1 — - and L2 (s) = — are strictly positive in ]0, C2[, it is assured that Bn (2, j) > 
for J = 1,2. 



Remark 3.3 As in Remark 3.2 h^Bn (2, j) when /i„ 0+. In particular, ft.„-B„ (2, 2) — > 
0. 

Analogously to the previous case, we present similar results for existence and uniqueness 
of nontrivial collocation solutions. 

Proposition 3.4 Let K he a kernel such that K{t,s) < K{t',s) for all < s < t < t' . 

Then there is existence near zero if and only if is unbounded. 

Proposition 3.5 Let K he a kernel such that K (t, s) < K {t' , s) for all < s < t < t' . If 

is unbounded near zero but it is bounded away from zero, then there is existence near 

zero. 

If, in addition, '^^^ is a strictly decreasing function, then there is at most one nontrivial 
collocation solution. 

For convolution kernels K{t, s) = k{t — s), the hypothesis on K means that k is increasing, 
and we can assure existence for fine meshes. 

Proposition 3.6 Let K he a kernel such that K (t, s) < K (t' , s) for all < s < t < t' . If 

is unbounded near zero and there exists a sequence {yn}n=i of positive real numbers and 

divergent to +00 such that lim„_j.-)_oo ^^^^^ = 0, then there is unconditional existence. 

If, in addition, is a strictly decreasing function, then there is at most one nontrivial 

collocation solution. 

As in the previous case, the nontrivial collocation solutions Zh are strictly positive (see 
Appendix A). Also, from the hypothesis on K and ([s]), it can be easily proved that the images 
under of the collocation points form a strictly increasing sequence, i.e. = Zq i < Zq 2 < 
Zi^i < Zi^2 < . ■ . < Zn^i^2- Particularly, if C2 = 1 then we can assure that z^ is a strictly 
increasing function. 



Remark 3.4 In the case C2 = 1, we can remove the hypothesis ^^K{t,s 
< s < t < t"' from Propositions 3.5 and 



3.6 



< K{t',s) for all 
because we do not need that Z12 > Zi i with 
I — I, . . . ,n— 1 for proving that z^ is strictly positive in ]0, t„[, and hence, we do not need to 
prove that Z„_2 > Zn^i (see Appendix A). On the other hand, if we remove this hypothesis, 
we can not assure that Zh is a strictly increasing function, as in the case m = 1 studied in 
Propositions |3.2| and |3.3| 



3.3 Nondivergent existence and uniqueness 

Given a kernel K, a nonlinearity G and some collocation parameters, we are interested in 
the study of existence of nontrivial collocation solutions using meshes Ih with arbitrarily 
small /i > 0. Following this criterion, we are not interested in collocation problems whose 
collocation solutions "escape" to +00 in a certain (7„ when /i„ — >■ 0"^, since this is a divergence 
symptom. 

Let S be an index set of all the nontrivial collocation solutions of a collocation problem 
with mesh I^. For any s e we denote by Zg.n,i the coefficients verifying equations Q (with 
n = 0, . . . , — 1 and i = 1, . . . , m) and such that, at least, one of the coefficients Zg-oA is 
different from zero (for some i € {!,..., to}). 

So, given K, G and some collocation parameters, we are going to define the concepts of 
nondivergent existence and nondivergent uniqueness of nontrivial collocation solutions. 
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Definition 3.1 Given = <o < • • ■ < such that there exist nontrivial collocation solutions 
using the mesh to < ■ ■ ■ < tn, we say that there is nondivergent existence in t^ if 

Zh^ inf <^ max {Zs-n.i} f 

seSh„ [_j=l,...,m J 

exists for small enough hn > and it does not diverge to +oo when /i„ — > 0+. 

Note that the index set of the nontrivial collocation solutions is denoted by Sh„ because 
if we change /i„ then we change the collocation problem. 

Definition 3.2 Given Ih = {0 ^ to < . . . < t^^i} such that there exist nontrivial collo- 
cation solutions using this mesh, we say that there is nondivergent existence if there is 
nondivergent existence in t+ for n ~ 0, . . . , N — 1. 

For studying the concept of nondivergent uniqueness we need to state the following defi- 
nitions. 



Definition 3.3 Given — to < ■ ■ ■ < tn such that there exist nontrivial collocation solutions 
using the mesh to < ■ ■ ■ < tn, we say that there is nondivergent uniqueness in tn if 

min <^ max {Zs-n.i} > = Zh„ 

s€Shn I 2 — l,...,m J 

exists for small enough hn > 0, and it does not diverge to +oo when /i„ — > but 

inf N max {Zs:n,i} \ - {Zh„} ] 

\ [_i=l,...,m J J 

diverges (note that inf ~ +oo). 

Definition 3.4 Given Ih = {0 — to < ■ ■ ■ < t^^i} such that there exist nontrivial collo- 
cation solutions using this mesh, we say that there is nondivergent uniqueness if there is 
nondivergent uniqueness in i+ for n = 0, . . . , iV — 1. 

When "nondivergent uniqueness" is assured, but there is not "uniqueness" (of nontrivial 
collocation solutions), there is only one nontrivial collocation solution that makes sense, as 
it is stated in the following definition. 



Definition 3.5 In case of nondivergent uniqueness, the nondivergent collocation solution is 
the one whose coefficients Z„ ^ satisfy max^^i^ „ {Z„ ,j} — Z^^ for n = 0, . . . , A'^ — 1. 

Next, we are going to study nondivergent existence and uniqueness for cases 1 (rn = 1 
with Ci > 0) and 2 (m — 2 with ci =0). Recall that the general conditions over K and G 
are always held, even if it is not explicitly mentioned. 

But first, we are going to state a result that reduces the study of nondivergent existence 
(and uniqu eness ) for any mesh Ih to the study of nondivergent existence (and uniqueness) in 
t^. Lemma 3.4 can be easily proved taking into account the general conditions over G. 



Lemma 3.4 Given a > 0, the minimum of the nonzero fixed points of G {a + j3y) exists for 
a small enough /3 > and converges to G (a) when /? — > 0+. 

//, in addition, '^''^^^•^ is strictly decreasing near zero, then the other nonzero fixed points 
(if they exist) diverge to -\-oo. 

As a consequence of this lemma, we obtain the following proposition. 
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Proposition 3.7 In cases 1 and 2 with existence of nontrivial collocation solutions, "non- 
divergent existence " and "nondivergent existence in tQ " are equivalent. 

If, in addition, '^^'^"^^^ is strictly decreasing near zero for all a > 0, then "nondivergent 
uniqueness" and "nondivergent uniqueness in tQ " are equivalent. 

Remark 3.5 The condition "^lii^-tM strictly decreasing near zero for all a > 0" is very weak, 

since lini„_j.o+ '^'•"^^•^ — +oo. For example, if G is twice differentiable a.e. without accumu- 
. y . . 

lation of non-differentiable points and without accumulation of sign changes of the second 
derivative, then this condition is held. So, if a collocation problem has nontrivial collocation 
solutions and, roughly speaking, G is "well-behaved" , then nondivergent uniqueness in t^ 
implies nondivergent uniqueness in any t^. 

Next, we are going to give a characterization of nondivergent existence, but first we need 
the following lemma, that can be easily proved taking into account the general conditions 
over G. 

Lemma 3.5 // is unbounded, then the minimum of the nonzero fixed points of G 
exists for a small enough /? > 0. In this case, this minimum does not diverge to +oo when 
/? — ^ 0^ if and only if ^^^^ is unbounded near zero. 

If, in addition, is strictly decreasing near zero, then the other nonzero fixed points 

(if they exist) diverge to +oo. 



Using Lemma |3.5| and Proposition 3.7 it can be proved the next result. 

Proposition 3.8 In cases 1 and 2 with existence of nontrivial collocation solutions, there is 
nondivergent existence if and only if ^^^^ is unbounded near zero. 

If, in addition, ^^"'^y^ is strictly decreasing near zero for all a > 0, then there is nondi- 
vergent uniqueness. 



Remark 3.6 If ^^^^ is unbounded near zero, then the condition "Mf|±l) strictly decreasing 
near zero for all a > 0" is very weak, using the same arguments as in Remark |3.5| So, 
if G is "well-behaved" (see Remark 3.5), nondivergent existence implies also nondivergent 
uniqueness. 



To sum up, combining Proposition 3.8 with Propositions 3.1 3.2 3.3 (case 1) and 3.4 



3.5 3.6| (case 2) and taking into account Remark 3.4 
• K {t, s) < K {t\ s) for all < s < t < t'; 



is unbounded near zero Nondivergent existence near zero. 



G(y) 
y 

Moreover, if G is "well-behaved" Nondivergent uniqueness near zero. 
• (Hypothesis only for case 2: C2 = 1, or K (t, s) < K {t' , s) for all < s < i < t') 



y is unbounded near zero but bounded away from zero Nondivergent existence 



G{y 
y 

near zero. 

Moreover, if G is "well-behaved" Nondivergent uniqueness near zero. 

For convolution kernels K(t, s) — k{t — s) ^ Nondivergent existence or uniqueness 
(resp.) for fine meshes. 

(Hypothesis only for case 2: C2 = 1, or (t, s) < K (t', s) for all < s < i < t'); 
^^Y^ is unbounded near zero and there exists a sequence {yn}nTi of positive real num- 
bers and divergent to +oo such that lim„_5.+oo ^^{p^ = Unconditional nondivergent 
existence. 

Moreover, if G is "well-behaved" Unconditional nondivergent uniqueness. 
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Figure 1: Case 1 (m = 1, ci > 0). Approximation yh (grey dots) for different stepsizes /i, 
with ci from 0.01 (lower) to 1 (upper). The solution y is also represented (black solid line). 



4 Examples 



4.1 Numerical study of convergence 

In this section we are going to show an example of a collocation problem (with nondivergent 
uniqueness in cases 1 and 2), and study numerically how the nontrivial collocation solution 
Zh converges to a solution of the implicitly linear HVIE (p|) when /i — 0+. In fact, we are 
going to show how the function yh = Vzh (see Remark 2.1) converges to a solution of the 
original nonlinear HVHIE ([T]). 

We are going to consider the equation 



v{t) 



{t-s){y{s)f'^ ds, ie[0,l]. 



(12) 



The kernel and the nonlinearity verify the general conditions and all the hypotheses of 

Hence there is unconditional existence 
I for assuring that the unique 



Propositions 3.3 (for case 1) and 3.6 (for case 2) 
and uniqueness; moreover, we can also apply Proposition [3 
nontrivial collocation solution is nondivergent. 



Since the unique nontrivial solution of the nonlinear HVHIE (12 1 is given by 



1 

i 

144 



we study the difference between this solution and yu when /i — > 0^ (see Figures [l]and[3| 



2! and ^ it is 

inerent stepsi 



represented the relative error 



ii\yh-y\ 



varying ci (case 1) or C2 



In Figures 

(case 2), for different stepsizes h. It is shown that the rate of convergence is the same as h. 
Moreover, a convenient choice of the collocation parameter can reduce the relative error more 
than two orders, as we see in Tables [T] and [2j Nevertheless, how to find a good collocation 
parameter is an open problem. 



4.2 Examples in other cases 

If the collocation problem is not in the scope of cases 1 and 2, the existence of nontrivial 
collocation solutions is not assured, even if the kernel and the nonlinearity satisfy all the 
mentioned conditions. 

For example, if we consider 

• K{t,s) = [t-s f , a > 0, 

.G(y)=^//^ &>1, 



10 




Figure 3: Case 2 (m — 2, ci =0). Approximation yh (grey dots) for different stcpsizes h, 
witli C2 from 0.01 (lower) to 1 (upper). The solution y is also represented (black solid line). 




0.3697 0.3697 0.3698 0.369: 



Figure 4: Case 2. Relative error varying C2 (from 0.01 to 1) for different stepsizes h. 
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Case 1 (to = 1, ci > 0) 


h ^ 0.1 


h = 0.01 


h = 0.001 


Max. rel. error 


2.1 

(ci = 1) 


3.2 • 10-1 
(ci = 1) 


4.1 • 10-2 
(ci = 1) 


Min. rcl. error 


2.5- 10-2 
(ci = 0.25) 


3.4- 10-3 
(ci = 0.175) 


3.4- 10-4 
(ci = 0.168) 



Table 1: Case 1. Maximum and minimum relative errors varying ci, for different stepsizes h 
(see Figure [2]) . 



Case 2 (to = 2, ci = 0) 


h = 0.1 


h = 0.01 


h = 0.001 


Max. rel. error 


4.2 • 10-1 
(C2 = 0.01) 


5.2 • 10-2 
(C2 = 0.01) 


5.4 • 10-3 
(C2 = 0.01) 


Min. rel. error 


1.5 • 10-3 
(c2 = 0.358) 


1.4- 10-^ 
(c2 = 0.369) 


1.4- 10-^ 
(C2 = 0.37) 



Table 2: Case 2. Maximum and minimum relative errors varying C2, for different stepsizes h 
(see Figure |4]). 



the general conditions over K and G are held, as well as all the conditions for assuring 
unconditional existence and uniqueness in cases 1 and 2. Moreover, by Proposition |3.8[ the 
unique nontrivial collocation solution is nondivergent. Nevertheless, in other cases we can 
not assure anything, as we will show in the next examples (where existence is considered at 
least in cto)- 

• Let 6 = 2. If TO = 2 and ci > then it can be proved that there is not any nontrivial 
collocation solution (for any a > 0). 

But, if m = 3, a = 2 and ci > then there exists a unique nontrivial collocation 
solution at least in ctq, independently of C2 and C3; on the other hand, if ci = then 
there is not any nontrivial collocation solution (contrary to what happens in the case 
TO = 2). 

Moreover, if to = 3 and a = 1 then there exists a unique nontrivial collocation solution 
at least in (Tqi independently of the collocation parameters. 

• Let 6 = 3, a = 1 and to > 2. It can be proved that there exists a unique nontrivial 
collocation solution at least in (Tq, independently of the collocation parameters. This 
is a special case, since the collocation solution for m > 2 coincides with the solution of 
the corresponding implicitly linear HVIE, z (t) — t/\/6. 



5 Discussion and comments 

A general theoretical analysis of numerical approximations of nontrivial collocation solutions 
for equation ([T]) is an outstanding problem which, to our knowledge, remains open. Among 
the reasons for such difficulty is the lack of uniqueness of solutions and the lack of Lipschitz- 
continuity of the nonlinearity, to name a few. 

Nevertheless, this work is a first step towards such general analysis: we give conditions 
for different kinds of existence and uniqueness of nontrivial collocation solutions aiming also 
for the convergence analysis (e.g. nondivergent existence). 

We have also studied numerically a concrete example showing the accuracy of the methods 



and how their errors depend on the collocation points (see Section 4.1 1. Therefore, in the case 
of nondivergent existence and uniqueness of collocation solutions, collocation methods have 
proved to be a valuable numerical tool for approximating solutions. Moreover, these methods 
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are strengthened by the fact that, if is a strictly decreasing function, the nonzero fixed 
points mentioned in Lemma |3.2| and Lemma |3.3| are attractors, so they can be easily found 
numerically with iteration techniques. Furthermore, if there is nondivergent existence, the 
minimum of such nonzero fixed points is also an attractor. 

It is worth pointing out that the conditions we have imposed over the nonlinearity in order 
to assure the nondivergent existence near zero of nontrivial collocation solutions, namely 
must be unbounded near zero, agrees with the conditions for existence of nontrivial solutions 
of the original equation (1). To be more precise, if is unbounded near zero then it is not 
Lipschitz-continuous, which is a necessary condition for the existence of nontrivial solutions 
of equation ([ij. 

A Proofs 

Proof of Lemma I3.1L 

(*)) Let us prove that if is unbounded (in ]0, +oo[), then there exists /?o > such 
that G {/3y) has nonzero fixed points for all < /3 < /3o- Let us take yo > 0, j/i := G (yo) > 
and /3q :— ^ > 0. So, given < /3 < /3o, we have 

G(/3yi) < G(/?oyi) = G(yo) =2/1, (13) 

since G is strictly increasing. Moreover, by hypothesis, there exists y2 > such that '^'^^ > 
^. Let us define j/s := ^, then 



G(/32/3) = G(j/2)> 1=2/3. (14) 



So, taking into account (131 and (14), G {/3y) has fixed points between yi and 7/3, since it is 
continuous. 

(-4= (i)) Let us prove the other implication. Given Ad > > 0, we take < /3 < < /So- 
Then, by hypothesis, there exists yo > such that G(/3yo) = J/o, ^nd so, taking yi :— f3yQ 

we have = i > M. 

yi p 

((i) {ii)) Let us prove that (i) implies (zi), for A > {). Let us take yo > 0, yi := 
G (^ + j/o) > and := ^ > 0. So, given < /? < and < a < ^ (we can suppose 
a > taking [3a < /3o), we have 

G (a + (3yi) <G[A + fiAyi) - G (A + yo) = 2/i, (15) 

since G is strictly increasing. Moreover 

G(a + /30) = G(a) > 0. (16) 



So, taking into account ( 15 ) and ( 16 ), G (a + /3y) has nonzero fixed points between and yi, 
since it is continuous. 
{{a) =^ (i)) Trivial. 



Proof of Lemma 3.2, Given a > and a small enough /3 > 0, it is easy to prove (taking 



analogous arguments as in the proof of Lemma 3.1) that there exists yi > such that 
G (a + /3yi) < yi, since is bounded away from zero. On the other hand, it is easy to 

prove that there exists ya > such that G {a + fSys) > ys, since ^^^^ is unbounded near 
zero (note that the case a > is trivial taking y^ — 0). So, G {a + (3y) has nonzero fixed 
points between y^ and yi, since it is continuous. Moreover, it is clear that if '^'^ is a strictly 
decreasing function, this fixed point is unique. 
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Proof of Lemma |3.3[ Analogous to the proof of Lemma [3. 2[ but taking any f3 > 0. 



Proof of Proposition |3.1| , 

(<J=) Let us prove that if '^^^ is unbounded, then there is existence near zero. So, we are 



going to prove by induction over n that there exist Hn > {n = 0, . . . , N ■ 
< hn < Hn then there exist sohitions of the system ([7| with Zq.i > 0: 



1) such that if 



For n = 0, taking into account Remark |3.2| and Lemma 3.1 (i), we choose a small 



enough Hq > such that < hoBo < Pq for all < ft-o < Hq. So, since the lag term is 
0, we can apply Lemma 3.1 -(i) to the equation Q, concluding that there exist strictly 
positive solutions for Zo,i. 

Let us suppose that, choosing one of those Zo,i, there exist Hi, ... , i?„_i > such that 
ii < hi < Hi (i = 1, . . . , n — 1) then there exist coefficients Zi^i, . . . , fulfilling 
the equation ([t]). Note that these coefficients are strictly positive by Remark 3.1 and 
hence, it is guaranteed that the corresponding collocation solution Zh is strictly positive 
in ]0,tn[. 

Finally, we are going to prove that there exists iJ„ > such that if < /i„ < iJ„ 
then there exists Zn^i > fulfilling the equation ^ with the previous coefficients 

Zq^I, . . . Zn~\X- 

Let us define A :— i^„(t„ + ci). So, taking into account Remark 3.2 and Lemma 
3.1/11^, we choose a small enough < H^ < 1 such that < h^Bn < /3a for all 
< /i„ < Hn. Then, we have < (^n,i) "Si Fn itn + ci) = A because the hypothesis 
over K (and the general eonditions), is strictly positive in ]0, and t„^i < t„ + Ci. 
Hence, we can apply Lemma 3.1(ii), obtaining the existence of Z„.i (that is strictly 
positive by Remark|3.1|). 



(=^) For proving the other condition, we use Lemma 3.1 -(i), taking into account Remark 



First we are going to consider a convolution kernel K{t, s) 



Proof of Proposition 3.2, 

k{t — s). Taking into account Remarks 2.2 3.2 and Lemma 3.2 we choose a small enough 
H > such that 

h [ \ ((ci - s) h) ds < /3o (17) 
Jo 

for aU 0<h<H. 

So, ii < h < H, we are going to prove by induction over n that there exist solutions of 
the system ([t]) with Zo,i > 0: 



For rt = 0, by ( 17) we have < hoBo < (3o, because < ho < h < H. Since the lag 
term is 0, we can apply Lemma 3.2 to the equation Q, concluding that there exist 
strictly positive solutions for .^0,1- 

Let us suppose that, choosing one of those Zq i, there exist coefficients Zi.i, . . . , Zn-i,i 
fulfilling the equation ([T]). Note that these coefficients are strictly positive by Remark 
|3.1[ and hence, it is guaranteed that the corresponding collocation solution z^ is strictly 
positive in ]0, tn[. 

Finally, we are going to prove that there exists Zn.i > fulfilling the equation ([t]) with 
the previous coefficients Zq^i, . . . Zn-iX- 

On one hand, taking into account ([6|, the lag term Fn {tn,i) is strictly positive because 
Zfi is strictly positive in ]0, i„[ and k satisfies the general conditions. On the other hand, 
< hnBn < f3o because < hn < h < H. Hence, we can apply Lemma |3.2| to the 



equation ain, obtaining the existence of Zn,i (that is strictly positive by Remark 3.11 
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If, in addition, ^^^^ is a strictly decreasing function, the uniqueness is an immediate 
consequence of Lemma [3?2| 

For general (nonconvolution) kernels, the proof is analogous but i?„ is given by ([8| and 
it does not only depend on /i„. The existence of /i„ lies in choosing a small enough Hn such 
that hnBn < /3o for all < /i„ < Hn- 



Proof of Proposition 3.3, The proof is analogous to the proof of Proposition 3.2 but in 



this case -ff > is any positive real number, and using Lemma 3.3 instead of Lemma 3.2 



Proof of Proposition 3.4 



(«^) Let us prove that if ^^^^ is unbounded, then there is existence near zero. So, we are 
going to prove by induction over n that there exist i/„ > (n = 0, . . . , — 1) such that if 
< hn < Hn then there exist solutions of the system ( 10 1 with Zq,2 > 0: 



For n = 0, taking into account Remark 3.3 and Lemma 3.1 (i), we choose a small 



enough Hq > such that < HqBq (2, 2) < /?o for all < ft-o ^ Hq. So, since the lag 
terms are and Zq i — G (0) — 0, we can apply Lemma 3.1-(i) to the equation (10), 



concluding that there exist strictly positive solutions for Zq 2- 

Let us suppose that, choosing one of those Zq^2, there exist Hi, . . . , Hn-i > such 
that ii < hi < Hi {i — 1, . . . ,n — 1) then there exist coefficients Zi,2, ■ ■ ■ , Zn^i,2 

Moreover, let us suppose that these coefficients satisfy 
1, and hence, it is guaranteed that the corresponding 



fulfilling the equation ( 10 ) 



Zi,2 > > for / = 1, . . . ,n 
collocation solution is strictly positive in ]0,t„[. 

Finally, we are going to prove that there exists iJ„ > such that if < /i„ < Hn then 



there exists Zn,2 > fulfilling the equation ( 10 ) with the previous coefficients, and 
Zn.,2 > Zn,\ > 0: 
Let us define 



A := Fn [tn +C2)+ k ((C2 - S)) Li (s) ds G {Fn (tn + Ci)) . 



So, taking into account Remark 3.3 and Lemma 3.1 Y**^ we choose a small enough 
< Hn < 1 such that < h^Bn (2, 2) < /3a for all < /i„ < Taking into account 
(|6|, the lag terms Fn {tn,i) are strictly positive for i = 1, 2, because Zh is strictly positive 
in ]0,t„[ and K satisfies the general conditions. Therefore, by ([9|, Zn.i — G {Fn (^n,i)) 
is strictly positive, because G is strictly positive in ]0, +oo[. Moreover, hnBn (2, 1) > 0, 
and hence, /i„-B„ (2, 1) Zn.i > 0. So, 



< [tn 



hnBn (2, 1) Zn,l < A, 



because K satisfies the hypothesis, Zh is strictly positive in ]0,t„[, the polynomial 
Li (s) is strictly positive in ]0,C2[, the nonlinearity G is a strictly increasing function, 

1,2. Hence, we can apply Lemma 



and in j — in 



for j 



(10 1, obtaining the existence of Z„_2- 



3.1 



(ii) to the equation 



Concluding, we have to check that Zn,2 > ^n,i- Since K satisfies the hypothesis in (j6|, 
Fn (in, 2) ^ Fn {tn,i), and hcncc, by the properties of G, we have 



G (f„ (t„,2) + hn E'=l Bn (2, j) Zn 



> G{Fn{tn,2)) > G{Fn{tn,l)) 
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For proving the other condition, we use Lemma 3.1 (i), taking into account Remark 



Proof of Proposition 3.5, First we are going to consider a convolution kernel K{t,s) = 
k{t — s). Taking into account Remarks 2.2 3.3 and Lemma 3.2 we choose a small enough 
H > such that 

(18) 



h I k{{c2 -s)h)—As< 
lo C2 



for all < h < H. From here, the proof is analogous to the proof of Proposition 3.4 where 
Hn = H and using Lemma |3 . 2| instead of Lemma 3.1 So, we do not need any A, and we do 
not need to check that F„ {tn.2) + hnBn (2, 1) Zn.i < A. 



If, in addition. 



G(y) 



is a strictly decreasing function, the uniqueness is an immediate 



consequence of Lemma |3.2 



For general (nonconvolution) kernels, the proof is analogous but i3„ (2, j) is given by ( 11 ) 
and it does not only depend on hji . The existence of hfi lies in choosing a small enough Hn 
such that hnBn (2, 2) < /3o for all < /i„ < i7„. 



Proof of Proposition 3.6, The proof is analogous to the proof of Proposition 3.5 but in 



this case 77 > is any positive real number, and using Lemma 3.3 instead of Lemma 3.2 
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